options = optimset('GradObj','on','Hessian','on');%梯度开启，二阶导数开启
[x,y]=fminunc('fun2',rand(1,2),options)
